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Abstract. To treat the spectral statistics of quantum maps and flows that arc 
fully chaotic classically, we use the rigorous Riemann-Siegel lookalike available for 
the spectral determinant of unitary time evolution operators F. Concentrating on 
dynamics without time reversal invariance we get the exact two-point correlator 
of the spectral density for finite dimension N of the matrix representative of F, 
as phenomenologically given by random matrix theory. In the limit N — ¥ oo the 
correlator of the Gaussian unitary ensemble is recovered. Previously conjectured 
cancelations of contributions of pseudo-orbits with periods beyond half the 
Heisenberg time are shown to be implied by the Riemann-Siegel lookalike. 
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1. Introduction 

Universal fluctuations in quantum energy spectra under conditions of full classical 
chaos are well understood in terms of Gutzwiller's semiclassical periodic-orbit theory 
[212]. The phenomenological description previously given by random- matrix theory 
has been fully recovered for individual (rather than ensemble averages of) systems 
from the unitary, orthogonal [215], and symplectic [3] symmetry classes. 

The analogous problem for periodically driven dynamics will be taken up in the 
present paper. A unitary Floquet operator F describing the evolution of such systems 
over one period of the driving and its quasi-energy spectrum will be investigated. 
Limiting ourselves to dynamics without time reversal invariance, we shall aim at 
recovering the spectral fluctuations of the circular unitary ensemble of random-matrix 
theory. Surprisingly, the semiclassical analysis will not only produce the two-point 
correlator of the density of quasi-energy levels in the limit of infinite dimension N of 
the Floquet matrix (which coincides with the correlator of energy levels of autonomous 
dynamics) , but even the finite-dimension variant (a periodic function of the spacing 
variable, with N times the mean level spacing as the period). 

Floquet operators turn out easier to analyze than Hamiltonians. It is the so 
called selfinversiveness of the secular polynomial of unitary matrices which comes 
to our help; it provides a variant of the secular polynomial, a "zeta" function 
£(<£>) = y / det(— F^)e~ lN(p / 2 det(l — e 1(p F), which is real for real values of the phase 
variable ip and enjoys a rigorous finite-iV Riemann-Siegel lookalike. Four zetas will 
be combined as Z = (C( 1 Pc)C( 1 Pd)C( 1 Pa)~ 1 C( ( Pb)~ 1 ) to a generating function Z which 
latter has the two-point correlator of interest as a suitable derivative. Due to the fixed 
order N of det(l — e lip F), the generating function is represented by finite Fourier series 
with respect to the variables (pc , <Pd the Fourier coefficients are related to the traces 
of powers of the Floquet operator, t n = TrF n , by the venerable relations given by 
Isaak Newton nearly 350 years ago. In contrast, the inverse zetas in Z have infinite 
Fourier series where the Fourier coefficients can again be expressed in terms of the 
traces t n . Inasmuch as those t n are determined by a Gutzwiller trace formula we 
get semiclassical access to Z\ the trace formula for t n is a well behaved one, without 
divergence necessitating regularization. 

When scrutinizing the semiclassical limit of Z according to the strategy just 
sketched we shall meet with a remarkable property, long hoped for but hitherto 
elusive. It is well to highlight that property right away in intuitive semiclassical 
parlance, before pointing to the rigorous algebraic origin. The finite Fourier series 
representing Z with respect to the variables (fc^D have the welcome consequence, 
due to the Riemann-Siegel lookalike, that periodic orbits with periods only up to half 
the Heisenberg time determine the constituent zetas C(vc) and C(v?d)- On the other 
hand, orbits with periods up to infinity enter the infinite Fourier series of the inverse 
zetas in Z. That 'asymmetry' in the pairs ipa,<Pb and ipcj^D makes for a vexatious 
inconvenience in evaluating Z. In previous work on autonomous flows the analogous 
difficulty was fought with the help of imaginary parts of the energy variables large 
enough to cut off the contributions of orbits with periods beyond half the Heisenberg 
time. In our present case, no trick at all is needed. Orbits with periods larger than half 
the Heisenberg time make exactly zero contribution to C(<Pc)C( l Pd) in Z. The algebraic 
basis for the property under discussion is the rigorous Riemann-Siegel lookalike for the 
zeta function. That powerful identity forces the coefficients of the secular polynomial 
det(l - e'fF) = Y!m=o A ne mv with n > N (and thus the coefficients of e mv in the 
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finite Fourier series for with |n| > -y) to vanish — even though Gutzwiller's 
periodic-orbit sum for the traces t n contains no information of the finite dimension of 
the Floquet matrix. 

The rigorous Riemann-Siegel lookalike can be carried over to the treatment of the 
energy level fluctuations for autonomous flows. One has to resort to a stroboscopic 
description and work with the Floquet operator for a suitable strobe period which 
secures a one to one relation between N energy levels and eigenphases. Again, finite 
Fourier series arise for the generating function Z for the phases ipc , ^>d ■ For the 
unitary symmetry class, the diagonal approximation yields, in the limit N — > oo, 
the exact Z known from the Gaussian unitary ensemble. As regards corrections we 
show that the present approach reproduces, for N — > oo, the previous results for all 
Wigner-Dyson symmetry classes. 

The rest of the paper is organized as follows. In Sect. [5] we present the 
Riemann-Siegel lookalike for the zeta function of finite dimensional unitary matrices F. 
Newton's relations are employed to express zeta in terms of traces of powers of F. In 
Sect. [3] the generating function is defined and decomposed in additive pieces according 
to the Riemann-Siegel lookalike. In Sect. |4]we invoke the Gutzwiller trace formula 
for powers of Floquet operators and show the diagonal approximation for the unitary 
symmetry class to produce the exact generating function known from Random Matrix 
Theory. Some hints are given on why we are certain the vanishing of corrections can 
be shown in a generalization of previous work on autonomous flows. In Sect. [5] we 
argue that autonomous flows can be described stroboscopically in terms of unitary 
Floquet operators. In the final section[6]we summarize the progress made through the 
rigorous Riemann-Siegel lookalike, comment on related previous work, and speculate 
on future ramifications. 



2. Riemann-Siegel lookalike for unitary matrices 

Towards capturing spectral fluctuations we look at the secular determinant 

N 

A(<p) = det(l - eW) = ^ A n e mv (1) 

where the coefficients A n include Aq = 1 and = det(— F). The determinant 
A(ip) is defined in the whole complex (p-plane. We may think of the determinant as a 
discrete one-sided Fourier transform of the coefficients with the inverse 

A n = r^A(<p)e-^. (2) 
Jo 2tt 

The information stored in A(tp) is also encoded in the (trace of the) resolvent, 

t(w) = Tr - . (3) 

To decode, a geometric series for the resolvent is useful, and that series takes different 
forms in the upper and lower halves of the complex (p-plane. Defining the 'traces' (of 
powers of the Floquet operator), 

t n = TrF n , n = 0,±l,±2,... (4) 

we have 

tM _/ EZo^tn lm^>0 

W ~ 1-E~=i e~ mv C lm^<0 ' (0) 
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In both cases, imaginary parts of the phase may be seen as convergence ensurers. 
As above for the secular coefficients, we recover the traces t n from the resolvent by 
inverting the one-sided Fourier transform, t n = |^ t(ip)e~ mip . 

Inasmuch as the determinant of any matrix equals the exponentiated trace of its 
logarithm we can express the resolvent in terms of the determinant as 

d_ < rv 

dtp 

The foregoing relation leads us to different expressions for the secular determinant for 
the two possible signs of Im<^>, 



t(<p) = N + ig-lnA(<p) o A(tp) = exp(-i dtp' (t(<p') - N)^. (6) 



An exp 



We have fixed integration constants by invoking A(cp) — > 1 for Im^ — > +oo in the 
first case, and A(ip) — > AmQ iNv for Im9? — >• — oo in the second case. By expanding the 
right members of ([7]) in powers of e v we formally get infinite series. However, the 
traces t n , t* n conspire to secure A n = for n > N. Even on the real ip-axis the two 
expressions exist and coincide, provided tp does not equal any of the eigenphases of F. 
The traces t n and the secular coefficients A n are directly related as well. The 

Taylor expansion of gives A(ip) = J2v=o (J2Zi ■ The secular 

coefficients are obtained by employing the multinomial expansion for the powers (•) 

and collecting the coefficients of e mv with n = 0, 1, To formulate the result, we 

consider all partitions of the positive integer n into sums of smaller positive integers 
I with multiplicities vi, namely n — X);=o i ^ v i- ^ or t Qe 'vector' v we define two 
properties, L(v) = J2i=o l ana ^ ^(^) = Sz=o l Vl - ^ e ^ nen nave [IHS]), 

{n=L(v)} 

V 1>1 

It is remarkable that the expression (JSj) holds independently of the degree N of the 
secular polynomial. Nevertheless, if the traces t n are evaluated for fixed iV we get 
A„ = for n > N. The expression ([5]) can be regarded as the explicit solution 
of Isaak Newton's recursive relation (N — n)A n = Ylm=o ^mAi-m f° r the secular 
coefficients; for that reason we shall refer to it as the Newton relation. 

The unitarity of the Floquet operator entails what is known as the 'self- 
inversiveness' of the secular coefficients [B], 

A n = A* N _ n A N . (9) 

That property makes sure that the set of the secular coefficients contains exactly N 
independent real parameters (as which one may see the eigenphases of F). For N 
odd, the A n with n = 1,2,...,^-^- together with An fully represent the secular 
polynomial since the remaining coefficients are given by selfinversiveness according to 
([9]). Similarly, for N even, the A n with n = 1,2,...,^ carry the full information. 
As a further consequence of selfinversiveness, the following variant of the secular 
determinant, to be called 'zeta function', 

iNip/2 



Cte) = y/A^e-^A^) (10) 
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is real valued for real tp. For the case of even N which we shall stick to from here on, 
it is convenient to split the zeta function in three additive pieces, 



A* N e 



f-i 

-if f A n e mv 

n=0 



n=0 



= C-(v) + Co + C+(v)- (ii) 

The terms C- an d C+ contain Fourier components e wtp with, respectively, v = 
— 1, —2, and v = l,2,...,f while Co is the zeroth Fourier component, 
Co = \J A* n An_ = /An_A% = Co- The decomposition (JTTJ) of zeta can be seen as 



a rigorous Riemann-Siegel lookalike [51-110] . 

An alternative form for the zeta function is useful. It is obtained by importing 
the representation (J?) of A(<p) to the definition (ITU|) . 



/A^exp { - if ^ - E~ i ^} Im^ > 

/^exp | + if 93 - J2n =1 Hr— } Im ^ < 



(12) 



That representation of zeta by an exponentiated infinite Fourier sum allows us to write 
the reciprocal 1/C as 



= $ V^expj+if^ + E^i^} lm^>0 



(13) 



The innocent looking sign changes relative to (jl2|) have in fact drastic consequences. 
Clearly, the Fourier series expansion of 1/C cannot be a finite-order polynomial in e ±1,p 
but rather is an infinite series, 



1 



Im (/? > 
Im ip < 



(14) 



wherein the Fourier coefficients A„ differ from the secular coefficients A n only by the 
absence of the sign factor in ©, 



A n 



E 



n 



l v 'v t \ 



(15) 



3. Generating function 

In analogy to standard practice for energy spectra we introduce the multiplicative 
combination of four spectral determinants, 

The angular brackets denote two averages. One is over 27r-interval of a real center 
phase </> defined through 

<PA/C = <P+^-, <Pb/d=<P jj-- (17) 

A second average is, in principle, necessary over a small offset window [BJ. The 
generating function Z deserves interest since it yields the two-point function of the 
level density through two derivatives, see [BJ and | Appendix A[ 
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To make the average (• ) = (2tt)- 1 / ^(0 wel1 defined, we must assign 
nonvanishing imaginary parts of opposite sign to the phase arguments in the 
denominator on the r. h. s. of (|16[) . Therefore, we endow the offset variables ca/b with 
positive imaginary parts. On the other hand, the phase arguments of the numerator 
zetas are unrestricted by the definition of Z. 

To prepare for the semiclassical analysis of the generating function Z we represent 
the two inverse zetas in the definition (|16l) by the infinite series (|14[) but employ the 
Riemann-Siegel lookalike (fTT|) in the numerator. The generating function thus becomes 
a sum of nine terms, 

z= z v = (c(<PA)- 1 a<PB)- 1 &{<Pc)c j (<PD)) . (is) 

i,j=+,-,0 

The symmetries 

Z+-(e A , e B , ec, e D ) = Z- + (e A , e B , -en, -ec) (19) 

Z ++ (eA, en, ec, en) = Z — (ca, en, -en, -ec) 

Z+o(e A , e B , e c , e D ) = Z 0+ (e A , e B ,-e D , -e c ) 

Zo-(e A , e B ,e c , en) = Z- (e A , e B , -e D , -e c ) 

show that only five components, say, Z y,Z ,Z-q,Zo+,Zoq are independent. 

Explicit expression for the Zij will be written as need arises. For now, we just note 
the example 

OO 2 1 

Z_+ =e WA+e B -ec-e D ) ^ ^ A a A* b A c A* d (20) 

a,b=0c,d=0 
X ^e^ a ~ b+c ~ d ^) e jf(ae A +be B +cec+de D ) . 

herein, the center-phase average yields the restriction a + c = b + d such that only the 
offset phases remain as variables. The summations over the integers c, d go from to 

Y — 1, due to the definitions of (±. Clearly then, Z |_ is a polynomial of order y both 

with respect to e~ le °'^ N and er leD / N , without zero-order terms in either variable. It 
is now convenient to introduce projectors VSl and V® leaving the powers Q lve c/ N an( j 
e iue n /N ^^jj < v < q unchanged and killing all others. Under the protection of 
the product Vi.T'R the sums over the integers c, d can be extended to go from to 
infinity, like the sums over a, b. With the help of the expansion coefficients (|8I15[) we 
get Z |_ represented by 

/ L(a)+L(c) = L(b)+L(d) 

Z +=V C<pDlj;(e A +e B -ec-e D ) ^ ^ 

\ a,6,c,d 

-S- r/t \a n +c n /f* \b rl +d n f_1 \c n +dn , \ 

TT \(—\ ^ ' e ijj(a n eA+c n ec+b n e B +d n eo) \ 

L\l\nJ KnJ a n \b n \c n \d n \ IT 

n—l / 

Similarly, we define projectors V§ , Vq which leave unchanged the phase 
independent term in the respective zetas while killing all e we c/ N anc ] e ive D /N v _^ 
0. We then have Z 00 = V§V${ ...), Z Q+ = V§V°{. . . ), and Z_ = V C V${. . . ) 
with the parentheses the same as for Z |_ above. The analogous series for Z reads 
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(- L(a) + L(c)+L(d)=L(b)+N 

Z__=V°V%le^ eA+eB - eo+eD) A* N ( 22 ) 

111 1 (111) v tl p i§-(aneA+c„ea+b„e B -d n e D ) 



l\l\n) \n) ajb n lc n \d n l 



Tl=l 

where the projector leaves unchanged the powers w j^j 1 q < v < W an( ^ 

kills all others. 

There is a price to pay for the introduction of the infinite sums over the vectors c, d 
in (I21I22[) . To ensure convergence we must allow for infinitesimal positive imaginary 
parts of e c / D , as already done for &aib- 

4. Semiclassics for maps 



4--1- Trace formula 

We now foray into asymptotics by invoking the Gutzwiller trace formula. For the 
traces t n of Floquet maps of periodically driven systems we have [6] 

t n ~Y-== £ e lS /\ (23) 

^ V|det(Af-l)| 

a sum over period-n orbits. The contribution of each orbit involves a phase factor 
where the action S (defined to include the Maslov phase and measured in units of 
Planck's constant) appears as the phase. That S is the 'time dependent' action which 
generates the periodic orbit in n steps of the classical map. The prefactor involves 
the primitive period n and accounts for the instability of the orbit in terms of the so 
called monodromy matrix M. For the sake of simplicity, we shall treat periodically 
driven single-freedom systems. We thus have a single positive Lyapounov rate A to 
reckon with for each periodic orbit. For beauty, we refrain from putting subscripts 
distinguishing different period-n orbits on the action and the monodromy matrix. 

The number J\f(n) of period-n orbits is known to asymptotically grow 
exponentially with the period n |lltil3j . On average over a small window An of 
periods centered at n, such that 1 <C An <C n, one has 

(Af{n)) = c" A 7n . (24) 

The finitcncss of the number of period-n orbits makes the trace formula ([2^]) a well 
defined object not in need of any regularization, quite in contrast to the more widely 
known one for the density of energy levels of autonomous dynamics (see next section) . 

Long orbits dominate the semiclassical asymptotics, and therefore we have written 
the trace formula (|23|) without regard for r-fold traversals of orbits with primitive 
period no = n/r. Assuming ergodicity we encounter one and the same Lyapounov 
rate for all long periodic orbits. The eigenvalues e ± ™ Ar of the monodromy matrix give 
the stability prefactor | det(A/ — 1)| -1 / 2 ~ e -«-W 2 SU ch that the powers of the traces 
in the Fourier series (j2~T|) and (p2|) for Z |_ and Z take the form 

© (n) =e-^"+ b "+ c "+^(£e^) ( E e"* s ) , (25) 
(^) Q " +C " + %) b " = e -^(-^-+-^)(Ee*f " +C " +4 (Ee-f". (26) 
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4-2. Diagonal approximation, unitary class 

We now focus on the diagonal approximation for the unitary symmetry class: only 
those terms from the multinomial expansion of the powers of periodic-orbit sums 
in (|25l26p are kept in which each phase factor e lS / h pairs up with its complex 
conjugate such that all dependence on the actions cancels. We thus immediately 

get (tn) m (C) 1 ™ $mm>- For fixed to, any choice of m action-phase factors e lS ^ h 

. di ? s .... 

comes with to! possibilities to pick the matching inverses. With any of these choices 

the TO-fold summation over periodic orbits leads to the TO-th power of the diagonal 

approximation for the product t n t* n of just two traces, 

C (Crl = <W m! MX ag = <Wm! {n 2 e- nXT N{n)) m , (27) 

- diag 

a behavior reminiscent of Gaussian statistics. 
Using the result (|2T|) for Z A ^' we face 

Z & }f = V C V D {...) (28) 
with the building block 

oo a n -\-c n — b n -\-d n = m , ]Ac n +d n 



(■•■)= ^ A+eB - e °- 6D) II E E (e- A ^ ( n)) m 

n—1 m—Q a n ,b ni c n ,d n 



We here do two binomial sums and sum over m, 

oo 

(...) = e WA+eB-ec-e D ) exp^ e - nAT 7V(n)(e I ^ eA - e'^ ec ) (e'^ es -e 1 ^ 6 - ) . 

n=l 

In the summand of the sum over the orbit periods n within the foregoing 
expression, the factor M{n) fluctuates about the average given by the HOdA sum 
rule (j2~i]) . Inasmuch as the cofactor depending on the offset variables e^/nfc/D varies 
smoothly with n on the scale An of validity of the HOdA sum rule, we can now replace 
the fluctuating J\f(n) by its average, e~ nXr N(n) — > 1/n. By this 'selfaveraging' in 
the semiclassical limit, the generating function becomes smooth without any average 
beyond the one over the center phase </>. 

Finally, using J2T S ~7r = — ln(l ~ e la: ) we ge10 

(...)= C i(^ + e B -e c -e D ) ° i 6 . ( 30 ) 

Let us first throw a glance at Z__ as given by (|22l26l27p . The phase matching 
required by the diagonal approximation yields the restriction a n + c n + d n = b n which 
contradicts the one obtained from the center-phase average, J2 n n ( a n + c n + d n — b n ) — 

\ Strictly speaking, we are not allowed to use the exponential proliferation factor for the average 
number of period-n orbits for small periods n. However, the error thus committed is of order 
(no/N) 2 with no some period above which the HOdA sum rule is reliable. For a rough argument, 
assume the offsets of order unity and consider ^Z™^ ( e ''' e ' 4 — e ' ™ ec ) ( e ' ™ BB ~ e'w^j/n ~ 

(eA — &c)( e B — e D)^2^i—i ~0Z ^ (~Ff) 2 '< ^ ne argument is easily modified to allow for the period 
2nN of the offset variables. 
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N. We conclude that Z^ 5 vanishes. Due to the symmetry (|19p we have Z^ s = as 
well. The full generating function thus comes out as 

z dia s = + (z d _% g + 4| ag + z^ s ) + (Z diag + Z di ag + z d i ag ) . (31) 

The first four terms are obtained from (. . .) by applying respectively the projectors 
VqVq ,T>£T>b ^-pC-pD ( pCpD . the rema i n i ng ones follow from the symmetry 

conditions (fT9")l . To let these projectors do their job we expand as (1 — e"^ 60 " 1 " 6 - ') -1 = 
S^Lo e'~&( ec+eD \ After elementary calculations given in |Appendix B we obtain a 
polynomial of order ^ in e^^ ec and e^~^ eD . Amazingly, the final expression coincides 
with the sum of (. . .) and its Riemann-Siegel complement (obtained by the replacement 

ec, &d ~> -zd, -ec), 

Z diag _ e -(e A +e B -e c -e D ) « ^ 

1 _ e lv( e A+eB) ± _ g— (e C +e D ) v ' 

. 1 — e -k( e A-ec) 1 _ p -f(ee-eD) 
, e 3(eA+e B +e C +eD) 1 ' ' 



Indeed, the denominator (1 — c"^ ( - ec+eD ' ) ) cancels from the sum of the two terms in 
(|32|) such that the polynomial character mentioned arises. 

Most remarkably, we have arrived at the exact generating function for the circular 
unitary ensemble [14] for finite N . In particular, the periodicity in the offset phases 
with period oc N is a nice finite-size effect for Floquet operators of periodically driven 
dynamics. 

Of course, in the limit N — > oc we get the generating function describing 
correlation decay at finite offsets 

lim Z diag = c ^A+e B -e c -e D )/2 (gA + e D )(e B + e c ) , 33 , 
n^oo (eA + es)(ec + en) 

_ ^{e A +e B +e c +e D )/2 ( e A ~ e C )(e B ~ e D ) 

{eA + e B ){ec + en) ' 

as for the Gaussian unitary ensemble. 

The case of odd dimension N leads to the same result. The pertinent calculation 
differs in technical details not worth being spelled out here. 



4-3. Off-diagonal corrections 

Inasmuch as the foregoing diagonal approximation already gives the exact result for the 
unitary symmetry class, all off-diagonal contributions must cancel. There is nothing 
much to add to the previous work for flows of autonomous dynamics as far as the limit 
N — > oo is concerned. 

Some nontrivial modifications are needed for finite N and the periodicity of Z. 
While a full proof of the cancelation of off-diagonal corrections will be reserved for a 
separate paper, we would like to at least give some hints here. 

As previously, quadruples of non-ordered sets of periodic orbits (pseudo-orbits) 
come into play such that the pseudo-orbits within each quadruple have action 
differences of the order of Planck's constant (see the next section for some more 
details). The (Feynman type) diagrams characterizing those quadruples remain 
unchanged, and so is the ordering of the diagrams in terms of the number of vertices 
(V") and links (L) by increasing values of L — V. The rules for associating analytic 
expressions with diagrams are modified. Straightforward computer-based summations 
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yield order- by-order cancelation of the corrections with L — V = 1,2, ...,6 which 
altogether involve 7404 non-equivalent diagrams. 

5. Stroboscopic maps for autonomous flows 

The rigorous Riemann-Siegel lookalike for finite dimensional unitary matrices can 
also be made available for autonomous flows, by choosing a suitable stroboscopic 
description. We propose to show that the status of previous work benefits from that 
modification. In particular, the diagonal approximation for the unitary symmetry class 
will again lead us to the exact finite- N result (|32[) . without needing large imaginary 
parts for the offset variables eciD- The periodicity in the offset variables is in this 
case an artefact of the stroboscopic description; only the limit N — ¥ oo is of interest 
here. 

5.1. Classical strobe period 

We consider a small stretch of the spectrum of a time independent Hamiltonian 
containing N ^> 1 levels {E\ < E2 < . . . < -Ejv} and the pertinent time 
evolution operator over a certain strobe period r, to be called Floquet operator, 
F = Ylu=i e~ lEf * T / h \ti)(\(J-\- Choosing the strobe period such that the eigenphases 
fn = -E/j 7 "/^" °f F j us t fill t ne interval [0 } 2tt] once, namely r = 2ttTiI(En — E\), the 
increasing order of the eigenenergies carries over to the eigenphases. For large N, the 
spectral fluctuations of the Hermitian Hamiltonian and the unitary Floquet operator 
become identical, since the neighborship of ifi and ippj loses any significance. 

Three time scales are of importance for us, the Heisenberg time Tjj = Nt, the 
Ehrcnfcst time Tg = A -1 hi(S/h) with A the Lyapounov rate and S some classical 
action scale, and the strobe period r. These times are ordered as t < Tg -C T#. In 
fact, we may take the strobe period to be a classical time, independent of Planck's 
constant; in view of Th ~ h~* +1 , classical r means N ~ ftT^ +1 . On the other hand, 
classical r means an overall energy span En — E% = 0(h), and therefore the mean 
level density therein, p = N/(En — E±), needs no unfolding. 

5.2. Trace formula 

We adopt the Gutzwiller formula for the density of energy levels [6lll5j. 



Here the sum is over periodic orbits on the energy shell E and S° is the energy 
dependent action (sometimes called reduced action) which as a generating function 
generates periodic orbits for fixed energy E; finally, T denotes the primitive period. A 
Legendre transformation connects the energy dependent action and its time dependent 
counterpart as S = S° — ET. We again employ det \M — 1| ~ e AT . The traces of the 
N x N Floquet matrix t n = J2k=i e~ mrEk ^ h can be written 




(34) 




(35) 
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The step functions select the energy interval [E — AE/2,E + AE/2] filled by the 
spectrum, such that tAE/Ti = 2n. Semiclassically for n > we have, 

i r E+AE/2 

t n ._LV dE , Te -.E'nT/H-XT/2^S°{E>)/H + e ->S°(E')/Hy (36) 

2lTh ^— ' Je-AE/2 

The reference energy at the mid-point of the spectrum will henceforth be set to zero, 
E = 0. To check how strongly the phase of the exponentials within the integrand 
varies across that interval, we expand [S°(E') - E'nT]/h = [S°(0) + (T - nr)E' + 
^S"(0)(E') 2 + . . .] /%. The third term in that expansion defines a width oc \ffi larger 
than the span of the spectrum such that the term can be dropped. The second term 
of the expansion, on the other hand, varies by A(j> = 2tt{T — nr)/r across the spectral 
range. We can conclude that the term e~ lS / h is negligible for n > while e lS / h plays 
no role for n < 0. Inasmuch as we need the traces t n only for positive n we have, 



\^ r L e iS /h-\T/2 J_ f AE/2 dE l e i(T-nr)E'/h 



t AE 



-AE/2 



T V T 



The function sine a; = sin7ra;/7rx favors orbits whose periods T do not differ from m 
by more than a few r (see next subsection). The sharpness of the energy and the 
uncertainty of the orbit periods of the order r are the important differences of the 
present trace formula to the one valid for periodically driven systems, see Eq. (|23[) . 
As a consequence of that difference it will turn out necessary to introduce a certain 
regularization below. 



5.3. Diagonal approximation, unitary symmetry 

We begin with calculating the product of t n t* n in the diagonal approximation following 
from the periodic-orbit expansion (13"T|) . Denoting by To the smallest period above 
which the average number of periodic orbits with periods in [T, T + dT is given by the 
exponential proliferation law dTe XT /T we have, 



-XT ■ 2 f T 

Miir I n 



T „f(|)^(l-„). (38, 

This expression formally diverges at the upper limit due to slow decay of the sine 
function at large values of the argument. It is here that the regularization announced 
just above becomes necessary. To that end, suppose that in Q35|) we smoothen the 
cut-offs at E' = ±AE/2 such that transitions from to 1 and back take place through 
intervals of the order of the mean level spacing SE — AE/N. That smoothing 
cannot noticeably change traces of the Floquet matrix. However, in the periodic-orbit 
expansion the contributions with periods significantly larger than the Heisenberg time 
will be suppressed^. The integral (J38J) then becomes finite with the main contribution 



We can, e.g., replace the thetas in J35H by the error functions i 



1 + erf 



/ E'-E±AE\ 



The 



2 

periodic-orbit expansion of traces <l 36 D will then acquire the Gaussian factor e ' h , Tjj = 

2nh/8E, and the result (139 1) will be correct to within corrections of order of IniV which may be 
neglected compared with relevant values of n. 
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from the region where the argument of sine is close to zero such that T can be replaced 
by nr. After that the integration interval can be safely extended to [—00, 00] with the 
result 



f°° dT . 2 






L~ smc 1 




= n 









(39) 



The diagonal approximation for products of higher powers of traces t™(t*^) k can 
be expressed in terms of [in^nldiag m the same way as for periodically driven systems, 
and we recover Eq. (j27]) . 

Likewise, all subsequent steps of the diagonal approximation in Sect. 14.21 are 
unchanged such that the generating function (|32p is rederived. In the present case of 
autonomous flows, only the limit of small offsets, Ca/b/c/dI^ ^ yielding (l33l) . is of 
physical interest. Two significant advantages over the derivation in [TJE1E] are worth 
noting. First, we now need only infinitely small imaginary parts of the offset phases 
e c/D- Second, as long as N is finite, the finiteness of the Fourier series in the variables 
e c/D is rigorously preserved; while the Riemann-Siegel components of Z individually 
have non- vanishing contributions from pseudo-orbits with periods T > Tjj/2, their 
sum enjoys cancellation of these contributions. 



5.4- Quadruples of pseudo- orbits for Z 

We recall the Riemann-Siegel lookalike (fTTjl and the ensuing decomposition (fT5|) of the 
generating function. It is now convenient to write the component Z |_ in the form 

00 f ± 



x exp 



I hn | e m vA _ gin^c ^ _|_ hi ^ e -» ii PB _ c m <po ~j I ^ 



=1 

7C -r>D 



which is easily seen to be correct when using the projectors V_ , V_ in (|12[) . 

After invoking the trace formula (|37|) the sums over the traces in (|40|) are expressed 
in terms of the periodic orbits, 



V hl & n V „ y jS°/n-\T/2 y J_ e mv sinc ( £ - n ) . (41) 



^— ' n ' ^— ' nr \ r 

n— 1 n— 1 v 

Inasmuch as orbit periods of the order of the Heisenberg time Nt dominate we may 
replace the factor 3L with unity, accepting a relative error of the order -A* . The lower 
limit of the sum over n can then be shifted to —00. 
Next, we employ the discrete Fourier transform 

oc 

e mv smc{x - n) = e la ^> (42) 

n— — 00 

where {tp} is the (27r)-periodic sawtooth function equalling ip in the interval [— n, it]. 
For use further below we note right away that restricting the sum over n to the range 
|n| < y means putting a soft cutoff on the parameter x to \x\ < because for fixed 
n and varying x the function sinc(a; — n) is weakly localized near x — n. At any rate, 
the sum in (|4"Tj) can be approximated as 
00 

h± c m '-P ^ y c iS°/h-\T/2+ 1 { V }T/r / 43 n 

n— 1 
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and herein the truncation n < ^ would entail the soft cutoff — < 4r. 

— r 
Four such sums appear exponentiated in the expression (|40[) for We Taylor 

expand all four pertinent exponentials and so obtain a fourfold sum over pseudo-orbits. 

Neglecting orbit repetitions and again representing the four phases PA/B/C/D as m 

(HTJ) we find [IJ[2j[6] 



pC-pD c i(e A +eB-ec-e D )/2 ^_-^ya+VD 



A,B,C,T> 

exp 



A (T A +T B +T C + T d ) /2+ 1 - (S A + S c - S B - S D ) 



exp 



- (T A {^ + ^} + T c {<j> +^}- T B {ct> - ^} - T D {0 - ^}) 



Here .4, £? (C,P) are the pseudo-orbits associated with the denominator (numerator) 
of the generating function; vq and vd stand for the number of orbits in C,T>. 

Next we turn to the center-phase average over the interval — ir <<))<%. We 
write the sawtooth functions as {4> + 6} = <fi + S — 2tt8(tt — <f> — S) and {(f) — 6} = 
cf> — S — 2ird(—ir — <fi + S) , for 5 > 0. The 4> dependence outside the step functions in the 
exponent is given by i0AT/r where AT = T A + Tc ~ Tb — Tjj is the period mismatch 
between the pseudo-orbits A U C and B U T>. For partner pseudo-orbits the period 
and action mismatches are both of order h, and therefore i^AT/t can be dropped 
in the semiclassical limit. The remaining elementary but cumbersome integral over 
4> is periodic in IMSJ°IIL_ The projectors V^VR serve to truncate the respective 
Fourier expansions at — y; the cumulative periods Tq/d of the pseudo-orbits C,T> 
are thus restricted to be smaller than half the Heisenberg time. For autonomous 
flows the remaining periodicity in the offset variables is an artefact of the stroboscopic 
description. Of physical interest is only the limit e A /B/c/D/N 0. In that limit, 
the step functions in the exponent can be dropped since they make only for O(j^) 
corrections. The projectors V^V® lose their meaning and are to be removed, but the 
soft cutoff for the periods Tq/d a t half the Heisenberg time must be kept. The final 
result thus reached, 

T c ,T D <T H /2 

Z _ + ^ e i(e A +e B -e c -e D )/2 ^_ -^yc + v D g~ A(T A +T B +T C +T D )/2 

A,B,C,T> 

x e i(T A e A +T B e B +Tcec+Toe D )/T H e i(S A + S c -S B -S D )/h 

is identical with the one previously obtained in 1,2,6 . In quite analogous a manner 

we recover, in our present context, the previous result Z ~ 0. Equally ignorable 

is then the distinction of the zeroth Fourier component £o and we can simply write 
Z = Z-+ + Z+-. 

We would like to note that the reasoning in the present subsection is not restricted 
to the unitary symmetry class. 



6. Summary and discussion 

We employ the rigorous Riemann-Siegel lookalike available for unitary N x N matrices 
to treat spectral fluctuations both for periodically driven and autonomous dynamics. 

For periodically driven systems, a finite Gutzwiller trace formula arises for the 
traces of powers of the Floquet matrix, t n = TrF™. Only period- n orbits are 
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involved in t n . The Riemann-Siegel lookalike then represents the spectral determinant 
by pseudo-orbits with cumulative periods up to n = N/2, half the dimensionless 
Heisenberg time N. — Restricting ourselves to the unitary symmetry class (no time 
reversal invariance) we find the diagonal approximation to yield the exact generating 
function of the two-point correlator of the density of quasi-energies provided by the 
circular unitary ensemble of random matrices |14j . for finite matrix dimension N. That 
generating function reduces to the one for the Gaussian unitary ensemble in the limit 
N oo. All corrections, due to partnerships of periodic (pseudo-)orbits generated 
in close self-encounters, must therefore cancel. Using computer based summation of 
the contributing partner pseudo-orbits we ascertain order-by-order cancelation in the 
lowest six orders; the order is given by the difference L — V, with V the overall number 
of relevant self-encounters and L their overall number of encounter stretches. 

For autonomous flows, we define a suitable strobe period and the pertinent time 
evolution operator to get a Floquet matrix capturing a finite number N of energy levels 
as quasi-energies. The pertinent trace formula, derived from Gutzwiller's one for the 
density of energy levels, needs regularization. Again, for the unitary symmetry class 
the diagonal approximation yields the exact finite- N CUE generating function and, 
with TV — > oo, the exact GUE result. We recover the previously obtained pseudo-orbit 
expansions p3-[3] for higher-order corrections in the limit N — > oo, for all Wigner-Dyson 
symmetry classes. 

A certain charm and, in fact, progress over the previous treatment of autonomous 
flows can be seen in two facts. First, the polynomial character of the secular 
determinant det(l — e ltp F) carries over to the generating function; the latter retains 
its polynomial nature under the pseudo-orbit expansion, even when the restriction 
imposed by the Riemann-Siegel lookalike, n < N/2, is lifted. Second, only infinitesimal 
imaginary parts are needed for the quasi-energy variables. That progress is owed to 
the rigorous Riemann-Siegel lookalike. 

We would like to pay respect to related work. The explicit solution of Newton's 
relations giving the coefficients A n of the secular polynomial in terms of the traces 
t n has a long history. Some recent references are [1H3]. The rigorous finite- N 
Riemann-Siegel lookalike has been noted before [5 HTU1ITB] . Already more than a decade 
ago, the possibility of semiclassically accessing spectral fluctuations using Newton's 
relations and Riemann-Siegel for unitary matrices was suggested 6,16 . Reasons for 
expecting cancelation of contributions from pseudo-orbits with periods beyond half 
the Heisenberg time in correlators of spectral determinants were put forth in |17j . 
More recently, such cancelation has been demonstrated for spectral fluctuations in 
graphs [18]. 

Postponed to future work is a demonstration of the cancelation of corrections to all 
orders of the 'diagrammatic' expansion for Floquet matrices of the unitary symmetry 
class. An adaption of the strategy of [2] should not be overly difficult to come by, even 
for the orthogonal and symplectic symmetry classes. 
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Appendix A. Two-point correlator 

The quasi-energy level density can be expressed using a 27r-periodic delta function, 
Pi<P) = E^ =1 %- Vm) = (27r)- 1 Er=-o O Ei 1 e l ^-^)-W0 + with real <p. The 
factor e _ l"l° fattens the delta function to an ordinary function, though one of 
arbitrarily small width. Subtracting the contribution = ~p of n — we have 
the fluctuating part of the level density 

n=l 

We define the two-point correlator as ((Ap(cp) Ap(ip' '))) where the angular brackets 
demand a double average defined as follows. The variables 93, y>' cover a square of area 
(27r) 2 in a (p^'-plane. In that square we have N 2 pairs of eigenphases ip^^ip^, such 
that an area (^r) 2 can be associated with each eigenphase pair. The average aimed 
at must produce a smooth function of the phase difference (p — tp' . We introduce a 
real center phase <fi and an offset e as ip — <fi + and ipf = <f>— jr. Both new variables 
range in 27r-intervals, < (f> < 2ir and — it < 4? < 7T. For the double average we choose 
an area of the form of a narrow rectangle in the (0-^)-plane, with widths 27r for the 
center phase <f> and some Ae/N for the offset. That rectangle must contain many pairs 
tp^,ip u , and therefore we must require 4p ^> jfe. That latter condition is a mild one, 
allowing the window ^ to be much smaller than the average level spacing. 

Now, the average over the center phase (■) = (2n)^ 1 dcj)(-) can be done and 
yields, with the help of the Fourier decomposition (|A.1[) . 

R(e) = -^Re(^|t„|V 2 "^-" 0+ ). (A.2) 

n=l 

The smoothing over the small offset window remains to be done. 

It is easy to ascertain that R(e) equals the real part of the complex two-point 
correlator of the resolvent 

= ( A - 3 ) 

where without loss of generality the offset variable can be taken with an arbitrary 
positive imaginary part. Inserting the geometric series for the resolvent and doing the 
center-phase average we get 

c ( e )-^(E^i 2e2me/Ar >- ( A - 4 ) 

n— 1 

Comparison with (|A.2|) shows that indeed R(e) = ReC(e). The complex correlator, 
in turn, is obtained from the generating function as d eA d eB Z 



e A/B/C/D~ 



1 / v^°° „i2en/JVu 



Appendix B. Diagonal approximation: action of projectors 

For the purposes of the present appendix it is convenient to shorthand as 

I h\ N / 2 

a = e le ^ N , b = e le »/ N , c = e le ^, d = e le °/ N , S = cd, F = ^> . 

1 — ab 
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The building block (...) in (1501) can then be written 



^- FS (1-*) 

OC 

= F5- N ' 2 [1 - (ad + bc) + abS] s " 



(B.l) 



We now write the components of the generating function obtained by action of 
the projectors V-,Vq on (...), namely, 

^diag _ y,C<pD ^ _ 



5"\-(ad + bc) E S")+ab[ 6V 



(B.2) 



-—N/2 



,u=-N/2 



, „=-(JV/2-l) 



and 







a 



c 

Zo+ S = V °V D (...)= -F- 



(B.3) 



The symmetries (fT9| involve the substitutions (c,d,5) — > (l/d, 1/c, 1/5) and provide 
three more terms of the generating function, 



^diag = p 



N/2 . . iV/2 N/2-1 

E^-f~ + ^)E^ E *" 



z diag = 



7 diag 
J +0 



-Fbc. 



The last remaining term without a Riemann-Siegel partner is 



ZZ s = V^(...)=F(l + ab). 



Collecting all these terms we obtain 

N/2 

rjdiag _ p 



N/2-1 N/2-1 

S"-(ad + bc) 6 v + ab 

v=-N/2 v=-N/2 v=-(N/2-l) 



= F5- N ' 2 
This expression is equal to 



(ad + bc) — ; — h abd 



1-5 



1-5 



1-5 



^diag _ p 



5 -n/2 (i - ad ) (i - bc ) , jtv^ I 1 - q/c) (i - &Ag) 



(1-5) 

(■•■)+ (• ■ •)l(ec,eD)^(— en,- ec) 



(1-1/5) 



which proves Eq. ([32 
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